ESTIMATING COMPLEX EIGENVALUES OF NON-SELF ADJOINT 
SCHRODINGER OPERATORS VIA COMPLEX DILATIONS 
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Abstract. The phenomenon "hypo-coercivity," i.e., the increased rate of contraction for 
a semi-group upon adding a large skew-adjoint part to the generator, is considered for 
ID semigroups generated by the Schrodinger operators —d^ + x^ + i'yf{x) with a complex 
f^ ' potential. For / of the special form /(x) = 1/(1 + Ixl"^), it is shown using complex dilations 

^N , that the real part of eigenvalues of the operator are larger than a constant times I'll 2/ ("^+2) 



CN ■ 1. Introduction 

_H ■ The purpose of this note is to take up up the following problem suggested by C. Villani 

o ,! im in the proceedings of the last ICM, quoted here (with notation slightly changed): 

r-| ! Identify sufficient conditions on / : IR — ^ IR so that the real parts of the 

'^ \ eigenvalues of 

as on operator on L^ go to infinity as |7| -^ oo and estimate the rate of 
^ I divergence. 

CN ! The stated problem is a model for the phenomenon of "hypocoercivity," with the specific 

JO I form motivated by applications described in [JJ and outlined briefly below. 
fT^ ■ It turns out that a sufficient condition on / for divergence of the real parts of the eigen- 

t~>^ . values may be obtained using a general result due to Constantin, Kiselev, Ryzhik and 

O I Zlatos |[2i Theorem 1.4]. Specifically, we will show below that: 

Theorem 1.1. The real parts of the eigenvalues ofLy diverge to cx) as j ^ oo if {x : f{x) = t} is 
essentially nowhere dense for each f G R. 

rS I Remarks. 1) A set S is essentially nowhere dense if S = S' UN where S' is no- where dense 

cd \ and N has Lebesgue measure zero. 2) Since the operators L^ are unbounded, we should 

specify their domains. Throughout we will consider bounded /, so the natural choice is 
to consider Lj on the domain of self-adjoin tness for the real part Lq: 



V = Ixp: I \d^\p\^ + x^\^p\^ < oo I 



It is then not hard to show that L-y is the generator of a contractive semigroup and that L;r^ 
is compact. Thus the spectrum cr[Ly) consists only of eigenvalues. 2) In fact the "if" of 
the theorem could be replaced by an "if an only if" ||3]] - a proof of the reverse implication 
would follow, for example, from the method used recently to derive a similar result in the 
context of diffusion with drift on a compact manifold BH. 

However, the question of estimating the rate of divergence seems to be much more 
delicate. In particular, it is likely that the rate will depend on the specific form of / and 
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thus we may not expect such a general result. Nonetheless, as we will show below, a 
number of interesting examples with analytic / may be analyzed with the technique of 
complex dilations, yielding, in one case at least, a very good estimate on the rate. 

The method of complex dilations, as applied below, essentially require that / be ana- 
lytic, possibly with a branch point. Furthermore, as the proof will show, the asymptotic 
growth of eigenvalues depends essentially on / only in the neighborhood of a critical 
point. If the potential has multiple critical points one may expect each of these to con- 
tribute. For the sake of clarity, we will not try to formulate a very general result, but rather 
focus on a function / with a single critical point to see how the behavior of the spectrum 
depends on the behavior of / near the critical point. Specifically, we will consider 

(1-1) fix) = ^ 

for arbitrary k > 0. The main result of this paper is 



Theorem 1.2. Fix k > and let f{x) be as in (|1.1|) . Then there is a positive constant Ck such that 
all eigenvalues ofLj lie in the half plane 

{z : Rez>CK|7l'/^'+'^}- 

Remarks. 1) Since the spectrum cr(Ly) consists only of eigenvalues, the theorem is suc- 
cinctly stated as the inclusion c^(L-y) C < z : Re(z) > C\j\'^^^'^~^-^^ L 2) The proof relies on a 
quadratic form estimate which only gives lower bounds. However, in |[1]] it is mentioned 
that numerical calculations by Thiery Gallay for k = 2 suggest that the obtained rate 7^^^ 
is sharp. This is to be contrasted with the rate j^^^ provided by the commutator methods 
ofVillanillJ. 

In fT], this problem is motivated as a simplified version of a spectral problem that arises 
in the linear stability analysis of a 2D hydrodynamic equation. In that context, one is 
mainly interested in understanding the t ^ 00 asymptotic behavior of solutions to the the 
initial value problem 

(1.2) dtxj^ix, t) = dlxp{x, t) - x^^p{x, t) - i^-—-—^p{x, t), ^p{x, t) = ^Pq{x), 

-L I A- 

The semi-group St = e~^^'^ generates solutions to (|1.2[) via tp{x,t) = Sttpo{x). It follows 
that the large t asymptotic behavior is governed by the eigenvalue of L^ with smallest real 
part. Indeed, 



e-f^ 



_ '^7,a 6 



for any a. < inf Re((r(L^)), where f or A C C we let Re(A) = {Re(z) : z G A}. Hence, the 
conclusion of the theorem implies that 



(1.3) / \it{x,t)\^dx < C^e" 



.C^2/(K+2)f 



for solutions to (|1.2|) . 

The estimate (|1.3|) is quite striking as, in some sense, the convergence of i/? — ^ is driven 
entirely by the "dissipative" term d^ip — x^xp on the right. Indeed, without the dissipative 
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term, the solution would be 



ip{xj) 



-iTt 



Wtpo{x) 



which has constant L^ norm. If the L^ conserving term — 17/(1 + |x|'^)i/7 is absent then the 
L^ norm goes to zero, but at the modest rate e~' given by finding the smallest eigenvalue 
of Lq. However the combination of the "conservative" and "dissipative" terms results in 
much faster convergence to zero. This phenomenon has been dubbed "hypocoercivity" 
by Villani, in analogy with Hormander's notion of "hypoelipticity" 

From another point of view, this is not so surprising. After all, those of us who add 
milk to our tea know to stir the cup of tea after adding milk to speed up the mixing of 
milk and tea. More or less this is the same phenomenon as what is described in eq. 11.21 
After all, on it's own the milk would eventually diffuse through the cup of tea. However, 
we may greatly decrease the time to acheive diffusion by stirring, which essentially adds 
a convective term to the diffusion equation for the density of tea. This convection on 
it's own, in an imaginary liquid with no dissipation, would only set the milk and tea in 
ever circulating motion - it is "conservative"! Together the dissipative and conservative 
motions combine into a flow which mixes much faster. 

In recent years the mathematical analysis of hypocoercivity and related problems has 
been the subject of much research. Constantin, Kiselev, Ryzhik and Zlatos 13 have an- 
alyzed the phenomenon from an abstract functional analytic perspective summarized in 
the following theorem: 

Theorem 1.3 (Theorem 11.4 of [2]). Let Hg be a self-adjoint, positive, unbounded operator with 
discrete spectrum and let F be self-adjoint and satisfy 



\F^\ 



< 



C {tp, Hotp) and 



e'^^ip,Hoe'^^ip 



< 



B{t){^,Hotp) 



with constants C and B{t) independent of tp and B{t) G Lj^^{0,oo). Then for any j > the 
operator Hq + ijF generates a contractive semi-group and the following are equivalent 

• For any t > 

-<(Ho+i7F) 



(1.4) 



lim 



0. 



The operator F has no eigenvectors in Q{Hq) = form domain o/Hq. 



Using Thm. [L3l we may easily prove Thm. lLlI Note that any normalized eigenvector \p 
satisfies 

g-t(Ho+i7f)w, = g-meA 

with A the corresponding eigenvalue. It follows that all eigenvalues satisfy 



(1.5) 



ReA> --In 

- t 



-t{Ho+iyF) 



Thus (|1.4[) implies that lim^-j.ooinf Re(7(Ho + 7F) = cx). This result applies in the present 
context, with the operator Hq = —d^ + x^ and F = f{x). The form domain of Hq is the 
space 

(1.6) ni = {xp:dx^par\dxip eL^iW.)}. 
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-t{Ho+iyF) 



^ Oas 



An eigenvector of F is a function supported on a level set of /. Thus 

7 — > oo if and only if the level sets of / support no (non-trivial) T-Li eigenf unctions. Since 
Til C C(R) Theorem O follows. 

However, (|1.5|) and (|1.4|) do not provide quantitative information on the rate of diver- 
gence. To find this we must estimate, for finite 7, the location of the eigenvalues. Before 
proceeding, it is instructive to ask, in general, "how do we estimate the location of spec- 
trum?" When one is available, a variational principle relating eigenvalues to extrema of 
a quadratic form is one of the most effective tools. For example, consider the self-adjoint 
eigenvalue problem 

(1.7) -dlxp{x) + V{x)tp{x) = \tp{x), 

with real valued V that diverges to +00 as x — > ±00. The smallest A G IR such that a 
solution tp G i-^(lR) may be found is called the ground state eigenvalue and satisfies a 
variational principle 

(1.8) Ao = infM \dxHx)\^ + V{x)\xp{x)\^dx / \xp{x)\^dx = l\ . 

[^ J —00 J —00 J 

Furthermore a minimizing wave function ip satisfies the eigenvalue equation (|1.7|) . The 
min-max principle generalizes this to higher eigenvalues. 

It would be hard to overstate the utility of (|1.8|) . To find an upper bound on Aq we just 
compute the energy £{tp) = J^^ \dxtp{x)\^ + V{x)\ip{x)\^dx of any normalized function 

xp G L^(R). To find a good upper bound we "simply" make a suitably clever choice of 
ip. Finding a lower bound is often more involved, requiring a uniform estimate from 
below on £{xp), but even that is often possible, particularly in an asymptotic limit (see, 
e.g., Simon's paper on semiclassical analysis [5]). In way of contrast, an exact analysis of 
the eigenvalue problem (|1.7|) is possible only for a few explicitly solvable examples. 

However, for the operators L^ considered here — which are not self-adjoint — a direct 
variational approach provides little insight. It remains true that any eigenvalue A satisfies 

(1.9) ReA > milRe{xp,LyXp) : f \'ip{x)\^ = l\ , 

since Re A = Re {ip\, L^ip^) with xp\ the corresponding normalized eigenf unction. How- 
ever, the infimum on the r.h.s. is insensitive to the imaginary term i7/(l + |x|'^)'" in the 
operator, and thus is independent of 7. So all we learn is that 

ReA > milRe{xp,Loxp) : f \xp{x)\^ = l\ = 1, 

using the well known explicit diagonalization of Lq = — 9^ + x^. 

The cornerstone of the variational approach is a relation between numerical range and 
spectrum valid only for normal operatorsjj The numerical range of a linear operator L on 
a Hilbert space H refers to the set 

num(L) = { {xp, Lip) y^ \ ip eT-L and \\xp\\-^ = 1} . 



Recall that L is normal if it commutes with it's adjoint L^. Lack of self-adjointness is not really the 
problem here. The variational approach would work in principle if Lj were normal. 
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If L is normal then we have 

(1.10) cch((r(L)) = cl(num(L)), 

where cch denotes the "closed convex hull" and cl denotes topological closure. This may 
be seen using the spectral theorem, since it is elementary for multiplication operators on 
L^ spaces. In particular, one has the following variational principle for normal operators: 
extreme points of the closed numerical range are in the spectrum. 

For a general closed operator L we do not have (|1.10)) . However, since the point spectrum 
of an operator L clearly falls in num(L) we have 

(1.11) a{L) C cl(num(L)), 

whenever L has compact resolvent — as do the operators L^ considered hereo In fact, 
since num(L) is convex [6^ Theorem V.3.1], we have 

(1.12) cch((7(L)) C cl(num(L)) 

for L with compact resolvent. However, there does not appear to be a relation between 
spectrum and numerical range valid for all operators L with compact resolvent beyond 



the inclusion (|1.12|) . Furthermore one knows from examples that the spectrum may lie 
arbitrarily far from the extreme points of the numerical range. For example, the numerical 
range of a 2 x 2 matrix is an ellipse with foci at the eigenvalues, but the radii of the ellipse 
may be arbitrarily large. 

One quite natural idea is to exploit the invariance of c^(L) under the map L i— > T^^LT 
where T is bounded with bounded inverse. This map fixes the numerical range only if 
T is unitary, so we may gain something by this procedure. Indeed, if we were lucky 
enough to find T so that T~^LT were normal, we would have a variational principle for 
the eigenvalues of L using the quadratic form (^xp, T~^LTy Even if T~^LT is not normal, 
we may hope to use T to bring the extreme points of the numerical range closer to the 
spectrum of L. In finite dimension, this procedure will work, in principle, for any L. 
Indeed, if dim(H) < oo then 

(1.13) cch((7(L)) = fl cl(num(r-iLT)). 

Tegcin) 

In fact, using the Jordan canonical form one can show a bit more: for any e > there is 
Te G QC{T-i) such that cl(num(r~^Lrfr)) is contained in the e neighborhood of cch((r(L)). 
It is not clear if (|1.13|) holds for general L with compact resolvent when dim('H) = oo. 



Anyway the answer to this question may not be so relevant, as in practice it is rather 
difficult to produce an effective conjugating operator without already knowing the Jordan 
form of L. 

So how do we locate the spectrum of L^? Let us first consider an heuristic approach, 
that is very close in spirit to the rigorous method applied below. Let A be an eigenvalue, 
Lrytp = Alp, with eigenf unction ip. Suppose that ip has a holomorphic extension ip{z) for 
z ^ Cr := {xe^^ : X G R and \9\ < j6}. By rotating the contour on which we evaluate the 



Eq. (|l.ll|l also holds whenever L is bounded, but may fail for an unbounded closed operator. For exam- 
ple consider Mip = ip' on the space Hq(0, 1) of functions vanishing at and 1 but with L^ derivative. Then 
M is closed and num(M) = ilR, but a{M) = C since ran(M — zZ) is perpendicular to e^^ for any z G C 



JEFFREY H. SCHENKER 



eigenfunction equation we find that 

-e-^''dl<pix) + e^'x^cp{x) + i^ ]^ cp{x) = X<p{x), 

-L I fc, -^ 



where (p{x) = tp(e^^x) for some 6 < ^. (We have extended |x|'^ to the holomorphic function 
(z^)2, single valued on C^ with a branch singularity at provided fi is not too large.) 
Assuming the complex rotated function (p is square tntegrable, we discover that A is an 
eigenvalue of 

(1.14) Lf) = -e-2iea2 + e^iV + i7 ^-^—. 

In particular, applying (|1.9|) with L^ ' in place of Ly, we find Re A > inf (r(Re L^^ ^ ) where 
ReLf = ^[Lf +Lf '] = cos(2e) j-B^ + x^ + ^-^^^'^^^ '^''^ 



The operator in curly brackets is a Schrodinger operator with (real valued) potential 

T./ N 7 stn(K0) Ixl'^ 



cos(20) M _|_ e'^i^lxKI^' 
If 6 is sufficiently small then V{x) > x^ and near the origin 

rrr ^ 2 sin(K0) , ,^ 

Applying ideas from semi-classical analysis, as in |[5l, one concludes for large 7 that the 

(16) 
ground state of Re Ly is an approximate ground state of 

^^n\ f -^2 sin(K0) , ,^ 
cos(29)|-a? + 7^kr 

whose ground state eigenvalue is seen to be proportional to 72/(^+2) |^y scaling. 

There are two deficiencies with the above argument. First we do not know that the 
wave function is analytic. Second, even if ip were analytic, there would be no a priori 
reason to believe that cp G L^. A reader familiar with the technique of "complex scaling" 
from scattering theory for (self-adjoint) Schrodinger operators will recognize the way out. 
We simply ignore these problems, focusing our attention instead on the spectral analysis 
of the following analytic family of operators 

(1.15) L^"^ = -e-2-a2+e2«'x2 + i^ L^ 

with w = u + 16 a complex parameter in a strip {|^| < 3} with S sufficiently small. For 
real w = u, with = 0, we have 

r (") T r Tt 

where Tu is the unitary dilation operator, 

(1.16) Tuipix) = e^Xe"^). 
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However, for complex w, Lr^ is not obtained from L^ by conjugation, but rather by ana- 
lytic continuation from the operator valued map u h^ Uy . With this set-up, we will prove 



Theorem II .21 in two steps: 

From analyticity w 

varies in the strip {| Imzf | < ^}. 

Using semi-classical analysis of 

mate on the real part of the lowest eigenvalue. 



(1) From analyticity we will show that the spectrum of L)^' is independent of w as w 

lid') 

(2) Using semi-classical analysis of Re L)y for 6^0 we will obtain an effective esti- 



2. Proof 

For each a > OletSa = {wgC : |Imw|<a:}. The first step of the argument is to 

1 
is a compact operator valued analytic map on a strip Sa^. 



show that w i— > 
Lemma 2.1. Let 



1—iry 



1 1 

a^ = TT mm . ^ - o 



IV) 



For each w G S^^ the operator L\ ' has compact resolvent and 

Proof. We will analyze Lj through the associated quadratic forms 

(2.1) q(!'')(V^) =e-2-||a,V^||J, + .2-11^^11 J, +i^ / U(x)\^dx, 

J — oo L -\- e X 

defined on the form domain "Hi (see (|1.6|) ). Our goal is to show, in the terminology of Kato 
lUl, that w H- )► Q)y^' is a holomorphic family of type (a), which is to say 

(1) w ^^ Q)y"'{^) is a holomorphic map for each \p G Hi, and 

(2) for fixed w the form Qx^ {■) is sectorial and closed on T-Li. 

Once we have shown (1) and (2), it follows by definition that Uy'' is a "holomorphic family 

of type (B)" and thus, by Theorem 4.3 of f6] that Ir^"' either has compact resolvent for all 

w or for no w. Since L^ = L!y has compact resolvent it follows that Uy" has compact 
resolvent for all w. 

Take (1) and (2) for granted for the moment. Then w i— )► L^y" is a holomorphic family 
with compact resolvent. By [6, Theorem Vn.1.9], we have the further dichotomy that 

A G C is an eigenvalue of Lx^ either for zy in a discrete set or for all w G Sa^. We may rule 
out the possibility of a discrete set by noting that 

Au+iQ) _ „tT(i0)7- 

L,ry ^U j Ul 

where Tu is the unitary dilation operator ((T]|. That is, the family u ^^ i^[u+w) [^ related by 

unitary conjugations and is thus isospectral. Thus Ly is isospectral as claimed and the 
lemma is proved once we verify (1) and (2). 
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Turning now to (1) and (2), note that (1) follows easily from the explicit expression for 

Qy . (Since Re e'^^" > on S^^ the integrand in the third term on the right hand side of 
(|2.1|) cannot develop a singularity.) To prove (2), namely sectoriality, note that 

where the dilation Tu clearly maps Hi onto Tii. Thus Q^' ^ is closed and sectorial if and 

(i9) 

only if Qy is closed and sectorial. 

do] 

To show that Qy is closed and sectorial, note that 

ReQ^^'\,p)=cos{2e){\\d,xp\\l2+\\xxp\\l2]-j r Im— -1^— |V^(x)|'dx. 

Thus, since cos{k6) > for \6\ < ix-k, we have 

(2.2) cos(2e)Qo(t/^) - |7| HV'IIl^ < ReQ^ ^(V^) < cos(20)Qo(V') + ItI MI2 > 

where 

Qo{ip) = \\d^ip\\l2+\\xxp\\l2. 
Now 

ImQ^'\ip) = sm{29) {- ||a,i/;||J, + \\xxp\\lA + ^ T Re J^ \xp{x)\^dx. 

>- -* J —CO i ~r c |X| 

Thus,by(|Z2l), 



ImQf)(i/.)| <sin(20)Qo(</^) + |7|||V'fL2<tan(20)ReQf)(V^) + |7| {^ + ^£(26)) "'^"'^' 

(16) 
which is to say Qy is sectorial. Furthermore, since Qq is a closed form with domain Tii 

it follows from (|2.2|) that Re Qy and hence Qy is closed on Hi. □ 

2.1. Analysis of Re Ly . Since Ly and Ly are isospectral, if A is an eigenvalue of Ly then 

ReA > cos(2e)min(r(H^(0)), 

where Hy{6) is the real part of — vm^'r ■ That is, 

TT /^N -^? ? sin(K0) Ixl'^ 



cos(20) |l + ei'^^|x|'^|2' 



Note that Hy{9) = H^y(—9). Hence without loss of generality we may take 7 > 0. For 
the rest of the proof, let us fix some 6 G {0,(Xk). Theorem 11.21 will follow once we show 
that the ground state eigenvalue Ag = ^o{0) of Hy{6) satisfies an estimate 



(2.3) Ao > Cj^. 



(In the end one could try to optimize over the choice of 9, however this would affect only 
the proportionally constant not the rate of divergence of the eigenvalue.) 
The Schrodinger operator Hy{9) is of the form 

-dl + V{x), 
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where 



V{x) 



^2 + ^sin(K0) 



\x\ 



cos(20) |l + ei'^^|x|'^|2 
Thus V{x) has a global minimum at x = 0, with 

s\n{ke) 



V{x) 



1- 



ixr + x2 + o(ixr^ 



x\ 



0. 



cos(20) 

Eq. (|2.3|) follows if we can show that the ground state eigenvalue of Hy{6) is asymptotic 
to the ground state eigenvalue of the anharmonic oscillator 

(2.4) Ka = -dl + a\x\'' 



with a: = 7 



sin(fce) 
cos(20) ■ 



Indeed, scaling shows that the ground state of K^ satisfies Ao{Ka] 



Let be a fixed C^ "cut-off" function, with 

(1) < <t>{x) < 1 for all X, 

(2) 0(x) =lfor |x| < 1, 

(3) 0(x) =Ofor |x| > 2. 

Fix V > 0, to be chosen below, and let 

/(x) = (t>[a.^x), with a = 7 



sm{ke) 



cos(20) ■ 

Define J{x) so that J{x)^ + }{x)^ = 1, that is J{x) = ^/l-}{x)^. By the IMS localization 
formula (see |7] Theorem 3.2]), 

(2.5) H^ie) = jH^ie)J+ jH^ie)} - {djf - {d]f. 

We will estimate each of the four operators on the right hand side of (|2.5|) separately. To 
estimate }H{9)}, we simply drop the kinetic term —9^, to obtain 



JH^{e)J > inf 

\x\>ix^ 



X^ + IX- 



\x\ 



;L_^giK0|^|K|2 



?> 



inf 



y^/'^ + cc 



V 



(1 + y)^ 



f. 



To minimize g{^j) = y^^^^ + oc ,^^ -.2 over y > a. ^^ we must compare 



a.' 



gia-"") = a-2^ + oi^:r- — z^ - Ca^-^" 



1 + oi- 



-KV\2 



with the value of g at the stationary points, which solve 

^yl-\l + yf + l=y. 

There are two solutions to this equation, i/i ^ 1 and 1/2 ^ Ca^r^ as a — > oo. Evaluating g 
at these points we obtain 

1 _2_ 

Thus, provided we take v < ^^z we conclude that 



jHy{e)j > C0i2+-Kf 



72. 
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as a 



oo. 



To estimate fHry(6)J we compare to the anharmonic oscillator (|2.4[) : 

1 



JHJe)J-JK,J 



x^ + a.\x\^ 



M _^giK0|^|;c|2 



1 



r 



> 



-a sup 

\x\<2a-^' 



\X\ 



1 



1+ \x 



\K\2 



f - 



-Cal-2-/2 



as a — > OO. Putting this together with the estimate for }Hy{6)y and the estimate Ka > 

2 

Ca2+K^ we see that 



H^(0) = }KJ + }{H^{e)-K^)} + JH^{e)}-{d,}y-(d,} 



\-lKVr2 



:r2 






.21/ 



Coc^-2kvj2^^^2v 



since |9x/|, |9x/| < Coc^. To minimize the growth of the error term (in brackets), let us 
choose V so that the two terms are of equal magnitude, namely v = jq^Sc' yielding 



H^{9) > Cj^-Cj^ 



and completing the proof. 
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